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This  technical  note  describes  a  microwave  circuit  simulator 
implemented  in  the  technical  computing  language,  MATLAB.  The 
simulator  takes  an  input  circuit  with  a  known  admittance  matrix, 
performs  a  node  swap  that  puts  the  external  nodes  (ports)  at  the 
upper  left  hand  of  an  equivalent  admittance  matrix,  then  reduces  the 
matrix  (regardless  of  size)  to  include  only  the  external  nodes  using 
Householder's  method.  The  resulting  2x2  admittance  matrix 
provides  all  the  information  to  uniquely  define  the  circuit  from  its 
external  ports  and  facilitates  straightforward  calculation  of 
Z-parameters,  S-parameters,  etc.  (In  our  case,  the  reduced  matrix  is  2 
X  2;  however,  the  method  applies  to  multiple  port  devices  as  well.) 
Although  it  can  be  used  by  itself,  the  simulator  is  designed  with  a 
bent  toward  empirical  (circuit  based)  transistor  models  and  may  be 
incorporated  into  small  and  large  signal  transistor  models. 


1.  Introduction 


This  technical  note  describes  a  microwave  circuit  simulator  implemented 
in  the  technical  computing  language,  MATLAB.  The  simulator  takes  an  in¬ 
put  circuit  with  a  known  admittance  matrix,  performs  a  node  swap  that 
puts  the  external  nodes  (ports)  at  the  upper  left  hand  of  an  equivalent  ad¬ 
mittance  matrix,  then  reduces  the  matrix  (regardless  of  size)  to  include  only 
the  external  nodes  using  Householder's  method.  The  resulting  2x2  admit¬ 
tance  matrix  provides  all  the  information  to  uniquely  define  the  circuit  from 
its  external  ports  and  facilitates  straightforward  calculation  of  Z-parameters, 
S-parameters,  etc.  (In  our  case,  the  reduced  matrix  is  2  x  2;  however,  the 
method  applies  to  multiple  port  devices  as  well.)  Although  it  can  be  used 
by  itself,  the  simulator  is  designed  with  a  bent  toward  empirical  (circuit 
based)  transistor  models  and  may  be  incorporated  into  small  and  large  signal 
transistor  models. 
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2.  Algorithms 


2.1  Simulator 

1.  Design  input  (circuit) 

2.  Assign  node  numbers 

3.  Determine  self  admittances  (diagonals  of  [Y]) 

4.  Deterinine  trans  admittances  (off-diagonal  elements  of  [Y]) 

5.  Frequency  plcm 

6.  Calculate  [Y(f)]  for  each  frequency  point 

7.  Apply  node  switch  algorithm 

8.  Apply  Householder's  theorem  to  obtain  reduced  admittance  matrix 
[Yred] 

9.  Calculate  other  parameters  e.g.  [S] ,  [Z],  etc.  from  [Yj^j] 

10.  Repeat  for  every  frequency  point 

2.2  Node  Switch  Algorithm 

It  is  desirable  to  reduce,  when  possible,  the  amount  of  data  involved  in 
performing  a  circuit  simulation.  One  approach  is  to  reduce  the  admittance 
matrices  so  that  only  the  external  nodes  (ports)  are  described.  Before  this 
can  be  accomplished,  however,  it  is  necessary  to  reconfigure  the  matrix, 
placing  the  external  nodes  in  the  upper  left  comer: 

1.  Assume  that  K  and  L  are  the  nodes  to  be  switched. 

2.  Start  loop  for  rows. 

3.  Start  loop  for  columns. 

4.  If  row  K  (K  or  L)  AND  column  K  (K  or  L) ,  then  Y'row,coi  =  Yrow.coi- 

5.  If  row  K  (K  or  L)  AND  column  =  K  then  ')('row,coi  =  Y'row,L- 

6.  If  row  K  (K  or  L)  AND  column  =  L ,  then  =  Y'^ow^k- 

7.  If  row  =  K  AND  column  =  L ,  then  Y'^o^v^oi  =  Yl,k- 

8.  If  row  =  K  AND  column  =  K ,  then  Y'^w^coi  = 

9.  If  row  =  K  AND  column  K  (K  or  L),  then  Y'^ow^coi  =  Y'l,co1- 

10.  If  row  =  L  AND  colunrn  K  (K  or  L),  then  Y',ow,coi  =  Y'k,co1- 

11.  If  row  =  L  AND  column  =  K,  then  Y'row,coi  =  Y'k,l- 
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12.  If  row  =  L  AND  column  =  L,  then  Y'row,coi  =  Y'k,K' 

13.  Repeat  for  every  frequency  point. 

Example 


^1 

^11  ^12  ^13 

^1 

h 

Y21  Y22  Y23 

^2 

h 

^31  ^32^33 

^3 
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h 

^11 

^13  ^12 

^1 

h 

= 

^31 

^33  ^32 

^3 

h 

^21 

Y23Y22 

^3 

2.3  Application  of  Householder's  Method 

Assuming  that  a  node  swap  has  been  performed  in  order  to  place  the  exter¬ 
nal  nodes  at  the  upper  left  comer  of  the  total  admittance  matrix,  Yjotai,  we 
can  now  reduce  this  matrix,  regardless  of  initial  size,  to  a  considerably  smaller 
one  that  entails  only  the  external  nodes  (1  and  2  in  this  example  for  a  two- 
port  circuit.)  This  is  achieved  by  virtue  of  Householder's  method  so  that  a 
matrix  can  be  broken  into  submatrices  as  follows  [3]: 


\y,oMh 


y.,  y„y„y,.y 


11  M2  M3  ^  14  ^  15 
21  ^22  ^23  ^24  ^ 25 
31  ^32^33^34^35 
41  ^42  ^43  ^44  ^45 


^51  ^52  ^53  ^54  ^55 


^11  ^12 
Y21  Y22 

^31  ^32 
^41  ^42 
^51  ^52 


^13  ^14  ^15 
^23  ^24  ^25 

^33  ^34  ^35 
^43  ^44  ^45 
^53  ^54  ^ 55 


Y  Y  • 

^  ee  ^  ei 

y.  y.. 

^  le  ^  ii 


where  subscripts  "e"  and  "i"  denote  external  and  internal,  respectively. 
Likewise,  Ohm's  Law  in  the  matrix  form  can  be  reduced. 


1  = 


h 

'^ee 

Ve 

^ie 
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We  recall  that  an  individual  element  of  the  admittance  matrix,  Y^l  can  be 
calculated  by  shorting  all  nodes  K  or  L  and  applying  an  internal  current 
voltage,  Vj,  and  by  observing  the  current  However,  assuming  a  passive 
device,  aU  the  internal  currents  are  zero  and 

I  =  Y  V  +  Y  y  • 

0  =  Y  V  +  Y  y  • 

y  • =-  Y"^Y  y 

i  ^  n  ^  le^  e 


Y  -Y  Y~^ 
^  ee  ^  le^  le 


V 


e 


Y 


reduced 


=  Y  -Y  Y~^Y 

^  ee  ^  le^  u  ^  le 


3.  Results  and  Comparison  to  ANSOFT  Serenade 


Listing  of  Code 


%TSjiice.m 

swap  of  nodes  k,r  'e.g.  1r'2,:L=  5 
%  starting  w/koows  Y{n,m)  ^ . 


5^50;pts=51;  CONSTANTS  ' 

y«[.2-j*.039rj’'.0398  0;  j*  0398  1-j*  0105  4;  0 -1  1.2]; 

.Rs=^ :  Rd’^O  ;Lg»0;  Ls*=0;  l.d=0; 

[n^]i=«L<!e(Y){  k=ss2;L-3; 


,ym0d(23;FR>=-li^.f '  ''  ''' 

*irni  Vf /  '•/■''  '' 

C-  'v<v''''  ^  V  ,  ''  '..i'S?*:'  ''  ;{?V  ,  '"'  " 

tcb.mufe  '  ,;;;::;v  ^  :' 

oase(<(i»w««fc)«^’row^i:>V&((cotoiiip'^^^  y)raodP<Riw^«taia-'R^=^^ 

cjM»n(i»vs^k)<SWiow$^i))’d£  '  :>^' "'imsdKiWiCflluii^  '  ' 

>  ,  case^t(row>^k)<S^row-f»L)>&«6i«^^  :  yiBpdP(rdw,<»l¥^>ll>^^^ 

'case(Kiw^&colainn«fetX  ‘  :/  ^  yi3»4l>(«w,^uinn4?R^»^^  '  "'" 

oas«(rdwp*s^t&coiiiw^  ,'  F  ’  ^  yinodH(i»w,coliMt»iM'^>»^  ' ', 

-e3se(rdw-^&;^coteinnr^^  ^  5!i«odi^mw,<folMiiv^^  ' 

case(ro^sF<^‘^.CCca3,t»W^^  3rao<a»{row,oPtttdin^R)«^0<j^  ' 

<ase(mw-“L&<»lunvJ^*la  W  *.  yniodP(^,adiih'i^i|JR^^^  ». 

,  casV<fow---^&,w1i^^  ^  ,  jf^ 

-fe:  ’'■'  ':’■  r ''^•'  '■'  -fir'  ',.  ■:'  '  ‘ 


Ifef  FR  l:p^J  ''-i  '  ''■;,Ct  ■  %  f  ^  \o'' 

%redace'^-^trixtpjijaidi9'fe3»^^»^^s)t3si^H^«s<*!i^d^Tibefi^m  '  - 
nii~2ltttm^2;  :-VV>"  >1'  4^M  .  '  ^  ^ 

yee(^;,FR)**3r^P^’t|»n,Rn?ia:,  .  ^;'  '^4v\5 'f 

yd(«:,FtRj|F*^'6aP(l^  .'N-;  ^  '  '' 

yic(./JRHm»3idP(Tt«‘«'laii,  .  >*/ 

‘  y»i'(:^]|rR)»«yniodP(iTO+l^atm  '  ^v-  .  f  '  '  ^ 

ut  m <^<|90del 

'  ^Ji'(FR)«f^3?ycd(ia«**5«^2;5«:-^^l^  .  lit 

‘ '  5imGa(i4,FKy-  3^P@d(2^3^Vdcy?i^  /  zn^!g2.|:ii^-^yrcdett23FRy< 

.ati»3d(^l,FR>--^d{2;i,FRyd<^  p  '^:\^zj!nbdC2^l^>»yfcdCI»rjRR^^d 


Listing  of  Code  (cont'd) 


%%%%%%?/o%%  7,  Calculate  extrinsic,Z,  paratneters  by  adding  pamsiiics  back  to  zmod  %%?/o%%%Vo%%% 

Zmnd(l,l,HR>  7mod(l,l,f'R)  t  (Rg'Us)  ^  j.*ww{KR).’*(Lg  t  L.s); 

Zmod(l,2»FR)-zmod(1.2,FR)  +  Rs  -j.*w\vfFR).*Ls: 

Zmod^2,l,FR)=zmod(2,l,FR)+  Rs  j.*ww(FR).*Ls; 

Zmod{2^dTlrzinod(2.24=1l)  +  (Rd+Rs)  +  j  *ww(FR) .♦(Ld  +  Ls); 

ZmodP(  1 , 1 ,FRV  Zinod(l ,1,FR)yZo;  ZmodH  1,2,FR)  =Zmod(l,2,FR)//o; 

ZmQdP(2 J ,FR)«2:mQd(2,l 4'R)./Zo;  ZmodP<2,2,FR>=Zmod(2.24"R)./Zo; 


%%%%%%%%%yM%%%%%^<>%%Y^Wo%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%  8.  Convert  Zmod  to  S_mod  %%%%%%%%%%% 

%%%%%%%%%y^>%%%%%%%%%%%%%%%%%%%%%%%%%%v,,%%%%%%%%%%%%%Vo%%%%v^%%%%yo%%%%yo% 

deltHFRK2modP(  1 ,1  ,FR>+ 1)  ♦(ZiTiodP(2a,FR)+ 1 )  -  (2modP(  1  a,FR).*ZmodP(2, 1  .FR)); 

Smod(1 , 1  ,FR)  (  (Zmt>dP(1,l  ,FR)-l).*fZmodP(2,2,FR)f  1)  -  (ZniodP(1 .2,FR).*/modP(2, 1  .FR))  )./deltl f FR); 

Smod(l  ;2,FR)*2  .“ZmodPt  1 ,2.FR)ydeltl  (FR); 

Smo^lA ,FR)-2  *ZmodP(2, 1  ,FR)^'dch  1  (FR); 

Smod{2,2,FR)-((ZinodP(I.l,FR)^1).'^(ZmodP(2,2,FR)-l)-(ZmodP{UFR))*(Zn5ixlP{2,LFR)))7dcitl(FR); 

Smodn(FR)  Sniod(M,FR);  Sm«dl2(FR>  Smod(l,2,FR);  Smod2!(FR)-Smt>d(2,l,FR);  Smod22(FR)  Smod(2.2,FR); 
end 


Conclusion/Recommendations 


Improvements  that  could  be  made  include  a  text-based  input  interface  simi¬ 
lar  to  a  SPICE  net  list  or  even  a  graphical  user  input.  As  implemented,  it  can 
easily  be  modified  to  be  included  in  larger  MATLAB  programs,  e.g.,  to  simu¬ 
late  an  empirical  (circuit  based)  model  of  a  transistor  in  order  to  compare  it 
with  measured  data. 
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